Quantitative modeling of spin relaxation in quantum dots 



O 

(N 

O 

Q 

in 



J. P. Hansen, 1,2 S. A. S0rngard, 1 M. F0rre, 1 and E. Rasanen 2,3 

'Department of Physics and Technology, University of Bergen, N-5007 Bergen, Norway 
1 Nano science Center, Department of Physics, University of Jyvaskyla, P.O. Box 35, FI-40014 Jyvaskyla, Finland 
''Physics Department, Harvard University, 02138 Cambridge MA, USA 
(Dated: December 16, 2011) 

We use numerically exact diagonalization to calculate the spin-orbit and phonon-induced triplet- 
singlet relaxation rate in a two-electron quantum dot exposed to a tilted magnetic field. Our scheme 
includes a three-dimensional description of the quantum dot, the Rashba and the linear and cubic 
Dresselhaus spin-orbit coupling, the ellipticity of the quantum dot, and the full angular description 
of the magnetic field. We are able to find reasonable agreement with the experimental results of 
Meunier et al. [Phys. Rev. Lett. 98, 126601 (2007)] in terms of the singlet-triplet energy splitting 
and the spin relaxation rate, respectively. We analyze in detail the effects of the spin-orbit factors, 
magnetic-field angles, and the dimensionality, and discuss the origins of the remaining deviations 
from the experimental data. 
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I. INTRODUCTION 

Quantum dots (QDs) are well-known examples of 
confined and quantized systems in semiconductor het- 
erostructures. Among other potential applications in, 
e.g., quantum optics, QDs have been experimentally real- 
ized as controllable quantum bitsi On the road to quan- 
tum computing applications it is essential to understand 
the electronic and spin-dependent relaxation processes in 

The coupling between quantized states and lattice vi- 
brations in semiconductor structures has been examined 
already for decades4^ In QDs the spin relaxation has 
attracted attention both theoretically^— and experi- 
mentally,!^^ Recently, Meunier et alJ£- measured the 
triplet-singlet relaxation rate in a two-electron (quasi) 
two-dimensional (2D) QD. They also derived a semi- 
empirical model reproducing the general trends of the ex- 
perimental data. However, a good agreement was found 
using a spin-orbit (SO) coupling that depended on the 
magnetic field, and that was significantly smaller than 
the SO strength reported in earlier studies^ 

In this work we describe a two-electron QD using both 
2D and three-dimensional (3D) models, respectively, and 
include (i) the description of a tilted external magnetic 
field, (ii) the ellipticity of the QD, and (iii) both the 
Rashba and the linear and cubic Dresselhaus spin-orbit 
(SO) interaction. We point out that none of the previ- 
ous theoretical approaches have aimed at such a complete 
description, which is shown to be important when com- 
paring with experiment at the most detailed level. The 
two-electron states are calculated numerically exactly by 
diagonalizing the Hamiltonian, and the phonon-induced 
relaxation rate is obtained from Fermi's golden rule. We 
find a reasonable agreement with the experimental re- 
sults^ provided that the 3D description and a moderate 
ellipticity of the QD are applied. The remaining discrep- 
ancy between experiment and theory can be caused by 
a more complex dot geometry of the experiment, and/or 



it indicates the limit of the widely used effective mass 
approximation for electrons confined by a harmonic po- 
tential. 



II. MODEL AND METHOD 

Our two-electron QD exposed to an external magnetic 
field B is described by the Hamiltonian 



H = 



i=l,2 



[h(n) + h so (Ti)] 



47re r e [ri - r 2 |' 



(1) 



where r 2 = (xi,yi,Zi) with i = 1,2 are the coordinates 
of the two electrons. The single-electron Hamiltonian is 
given by a sum of the kinetic, external potential, and 
Zeeman terms as 



(2) 



+ -m* [Sujqx 2 + (l/S)Lu y 2 + u z z 2 } + <?> B B • S. 

We consider magnetic fields tilted from the xy plane with 
an angle 8, and tilted azimuthally from the y axis with 
an angle <f>. Thus the magnetic field has an expression 

B = Bo (cos <f> sin 8, sin <j> sin 8, cos 8). (3) 

The vector potential in the symmetric gauge, A = 
(B x r) /2, then becomes 

A = (z s'm 8 s'm (f) — y cos 8 , 

x cos 8 — z cos 4> sin 8, 
y cos <f> sin 8 — x sin <f> sin 8) , 



(4) 



which approaches A = Bo(—y, x)/2 in the 2D limit (z — > 
0). 

The harmonic confinement potential is asymmetric 
both laterally and vertically: the ellipticity in the xy 
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plane can be tuned with the (^-parameter, and the in- 
plane and off-plane confinement strengths are fixed to 
htjjQ = 2.5 mcV and haj z = 11.9 meV, respectively. The 
material parameters are chosen to be those of GaAs, 
i.e., effective mass m* = 0.067m e , relative permittivity 
e r = 12.4 and the gyromagnetic ratio g* = —0.44. The 
constant /ib = e/(2m e ) in Eq. <j3j) is the Bohr magneton. 

A harmonic form has been shown to describe well the 
electron confinement in both lateral and vertical semi- 
conductor (GaAs-type) quantum dots^ The model was 
established soon after the first Coulomb-blockade experi- 
ments, when the measured and calculated addition ener- 
gies were compared.— More recently, the harmonic model 
has been shown to explicitly yield the measured singlc- 
eletron spectrum^ and also the spin-blockade oscilla- 
tions for both lateral and vertical QD devices up to about 
50 electrons.—"^ 

The SO interaction in Eq. ([lj is given by 



Hr + hi 



with the Rashba term 



(5) 



(6) 



and the Dresselhaus term 2 ^ containing both linear and 
cubic contributions, 



hi 



7 



(vl){PyV y ~Px<J X ) 
l^(jPxP 2 yO X ~ PyPl° V ) 



(7) 



The coupling strengths of the Rashba and Dresselhaus 
terms are set by parameters a and 7, respectively, and 
a x and a v are the Pauli matrices. We choose a = 
and 7 = 11.15 eVA as our default values, but consider 
also different magnitudes to assess the sensitivity of the 
spin relaxation to the SO coupling. The measurements of 
Mcunier et aZ.— and 2D calculations^^ 3 - suggest relatively 
small coupling constants. This is in agreement with the 
experiment of Zumbiihl et al^ reporting 7 = 9 eVA 3 
which is notably smaller than the commonly used value 
of 27.5 cV A 3 . 

The total wave function is a product of the xy (pla- 
nar) component ^(ri,^), obtained from exact diag- 
onalization of the Hamiltonian in Eq. ([1]) in a (large) 
basis of spin-symmetrized products of Cartesian single- 
electron harmonic-oscillator eigenstates ; 22 i 26 and a sim- 
ple z (perpendicular) component ip z (zi,Z2) incorporat- 
ing the material thickness. The z component of the wave 
function is conveniently modeled in terms of the ground 
state wave function of a harmonic potential, with con- 
finement strength uj z defining an effective confinement 
length, L z = Ah/ (m*uj z ), in the z direction. 

The spectrum is obtained by diagonalizing the Hamil- 
tonian in a basis of spin symmetrized two-electron 
harmonic-oscillator states. Basis sizes \n x ,n y ,n z ) up to 
n-max = 6 for the (x, y) coordinates and n z = for the 



vertical coordinate are included. Analytical matrix el- 
ements are then obtained for all terms except for the 
electron-electron interaction in Eq. (JT]). 

It was shown by Popsueva et ali^ that, in the 2D case 
with lu x = uj y = wo, all matrix elements can be evalu- 
ated analytically. In particular, for any double pairs of 
basis functions the four-dimensional integrals over 
the electron-electron interaction can be expressed as 



1 2D 



,1/2 



where 



E 



(8) 



(9) 



The variable s M is here defined from the Bethe integral, 
1 1 f d 3 s 



l r i -r 2 



2tt 2 



(10) 



When this expression is introduced in calculating the ma- 
trix element between any pair of basis functions (i, j), the 

expression above is derived with s M = A /s 2 + s 2 When 



the integrals have been analytically evaluated once, the 
matrix elements for any confinement ujq are obtained by 
a simple multiplication. 

In 3D the electron-electron repulsion amounts to more 
complicated six-dimensional integrals. When extending 
with harmonic confinement also in the z direction a sim- 
ilar scaling relation can be obtained, i.e., 



T r 12 _ 
J 3D " 



1/2 



1-erf 



1 

2^, 



x e^zds,. 



(11) 



Here erf(x) = 2/y/n J x exp(— t 2 )dt is the error function. 
The matrix elements for the electron-electron interaction 
can thus be obtained for any ujq by scaling as in 2D when 
they have been calculated for a given value of the verti- 
cal confinement. We note also that the formula above 
provides a continuous route from 3D to 2D, since 



r 2 E = J™J r 3 n(^)- 



(12) 



For 2D calculations the finite thickness can be consis- 
tently taken into account by replacing Eq. © by Eq. 
(fTTj) when calculating the interaction terms. 

The triplet-singlet relaxation rate is calculated using 
Fermi's golden rule, which - in the case of phonon emis- 
sion - may be written asj2£ 



(5|F ph |T) 6(AE-hc ), 



(13) 
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where T and S refer to the initial singlet and final triplet 
states, respectively, V is the normalization volume, AE 
is the singlet-triplet energy splitting, q = (q M ,g z ) is 

the momentum of the released phonon, and i? p h = 
Si=i 2 e -l< - q "' ri+,3iZi ' is the phonon coupling. Linear dis- 
persion relations have been adopted, i.e., we set w§,t = 
AE/h = c ht q. Here q = 4720 m/s (c t = 3340 m/s) 
denotes the longitudinal (transverse) speed of sound.— 
The electrons couple to longitudinal acoustic phonons 
through the deformation potential couplin g 28 ' 29 



|Mi(q) 



2pc\V 



as well as the piezoelectric coupling, 

2 32 A 2 fo 2 4 (3q x q y q z ) 2 
e 2 r pciV |q| 7 



(14) 



(15) 



whereas coupling to transverse phonons only takes place 
through the latter one, 



32tt 2 he 2 h 2 14 (q 2 x gl 



q 2 y q 2 z + q 2 z q 2 x (3q x q y q z ) 2 



e 2 pc t V 



(16) 

In these expressions p = 5300 kg/m 3 is the GaAs mass 
density, Sa = 6.7 eV is the deformation potential con- 
stant, and /ii4 = 1.4 x 10 9 V/m is the piezoelectric 
constant. There are two transverse phonon modes, and 
hence Ms(q) is considered twice when computing the 
rate. 

The matrix element (S , |i? p h|T 1 ) in Eq. (jTB")) separates 
into a product of a planar (xy) and a perpendicular (z) 
component, i.e., 

( s ')| - i q„- r ,u/,( T )\ 



(S\H ph \T) = S) \e- i ^^r){^- iq ^ z ) 



i=l,2 



k,l.k' 



HI 



4m*w, 



Here q x = hq x /(2m*uj ) 1 / 2 , N = min(n,n / ), M = (|n - 
n'\ — T])/2, and r\ = 1 if n + n' is odd and otherwise 
(analogously in y)22. The phonon induced triplet-singlet 
relaxation rate is now calculated by inserting Eqs. (fT4"]) - 
(fT7|) into Eq. (fi~3|) and performing the integration over 
the phonon momentum q. 



III. RESULTS 
A. Singlet-triplet energy splitting 

Figure QJa) shows the total energies of singlet (thick 
lines) and triplet states (thin lines) in a two-electron, el- 
liptic QD as a function of the magnetic field. The clliptic- 
ity of the QD is fixed to 5 = 1.3 [see Eq. ©]. The triplet 
states with S z = 0, ±1 are split into three curves due to 
the Zeeman effect which is linear in B. We show results 
for three tilting angles of the magnetic field, 9 = 65° 
(outermost lines), 50° (middle lines), and 35° (innermost 
lines) as indicated in the figure. 

The tilting angle 6 has a relatively large effect on the 
energies, so that the closing point of the singlet-triplet 
energy gap AE moves to larger magnetic fields when the 
angle is increased. Figure [TJb) shows the gap as a func- 
tion of B in comparison with the experimental values 
(circles). First, the downward "bending" of AE at small 
fields is due to the ellipticity; in a circular QD the gap 
increases linearly as B approaches zero. In this respect, 
we find a good agreement with the experiment, and may 
conclude that the ellipticity of the real QD device is close 
to our chosen value 5 = 1.3. Secondly, the best overall 
agreement is found with a tilting angle 6 = 50° which is 
smaller than 9. 



(17) 



Here If/ 1 ,, ) and IV 1 ,, ) are the planar two-electron sin- 
glet and triplet states, respectively. We have assumed 
that the z component of the wave function is frozen, i.e., 
no excitations are considered in that direction. The ex- 
pansion coefficients akj are obtained from the diagonal- 
ization of the Hamiltonian in Eq. (TTJ) . njt and are the 
single-particle harmonic quantum numbers of one of the 
electrons in the x and y direction, respectively. The fac- 
tor of two in Eq. ((TTJ) arises from the fact that the total 
wave function of each state is antisymmetric. The single- 
particle matrix elements $ ni ,n./ = (4>ni \e~ lqxX \4 > n i ,) are 
essentially associated Laguerrc polynomials in q 2 multi- 
plied by a damping exponential function, 



(-ir^f x ^e^(-ql/2)L% M ^(f x ). (18) 
yn\n\ 



exp 



68° ± 5° reported in Ref. [If 



Another angle affecting the energy gap is the azimuthal 
direction of the magnetic field with respect to the axis of 
the ellipticity. In Figure He) we plot the variation of 
AE(B) as a function of the azimuthal angle (j) for three 
magnetic field strengths. We observe that the variation 
increases with increasing magnetic field, so that close to 
the crossing (B ~ 3T) the variation is ~ 20%. Combin- 
ing this observation with Fig. [TJb) leads to a conclusion 
that the best agreement with the experimental result is 
obtained with </> = 90°, i.e., when the shortest axis of 
the ellipse is parallel to the tilting direction of the mag- 
netic field. In the following we set <j> = 90° unless stated 
otherwise. 

With a fixed 5 obtained from the bending as B — >• 
[see Fig. Hfb)], the singlet-triplet energy gap essentially 
depends on loq and uj z . SO parameters are so small that 
they play only a minor role in the energy gap. The pa- 
rameter values used here, Hluq = 2.5 meV and fku z = 11.9 
meV, give a reasonable agreement with the experiment in 
the energy separation and in the spin-relaxation rate (see 
below). An increase in both parameters and thus smaller 
vertical extension can also lead to good agreement in the 
energy gap, but a poor agreement in the spin relaxation 
rate. This effect is further illuminated in Sec. III.C. 
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FIG. 1: (Color online) (a) Energy of singlet ground state 
(thick lines) and three lowest triplet states (thin lines) in a 
two-electron harmonic (Run = 2.5 meV, hw z = 11.9 meV), el- 
liptic (8 = 1.3) quantum dot as function of the magnetic field 
strength for three different tilting angles of the magnetic field, 
9 = 65° , 50° , and 35° . In the triplet states the dashed, solid, 
and dash-dotted lines indicate S z = —1, 0, and +1, respec- 
tively, (b) Singlet-triplet energy splitting for the three cases of 
the upper panel. Open circles correspond to the experimental 
data in Ref. [l5|. (c) Variation of AE(B) as a function of cj> 
for B = 0.5 T (solid line), 1.5 T (dash-dotted line), and 2.5 T 
(dashed line), respectively. The results for B — 1.5 T and 2.5 
T have been multiplied by factors 1.4 and 2.4, respectively, 
to set the scale. 



Wc point out that the results are not sensitive to us- 
ing other types of frozen ground-state vertical basis func- 
tions and/or a larger number of vertical basis functions. 
The latter, as well as representing the vertical dimen- 
sion by an infinite barrier at z = and a harmonic con- 
finement for z > does not lead to significant changes. 
We conclude that good agreement with the experiment 
is obtained for a narrow range of harmonic confinement 
parameters. However, the theoretical magnetic tilting 
angle needs to be reduced to 9 = 50° to arrive to this 
level of agreement. The origin of this discrepancy pose 
a challenge for future experiments and theory. In the 
remaining part of this work we will stick to 8 = 50°. 



B. Spin relaxation rate 

In Fig. [5] we show the calculated total relaxation rate 
(solid line) in comparison with the experimental result 
(circles with error bars) in Ref. UH. Dashed and dash- 
dotted lines correspond to the contributions from the 
piezoelectric coupling and the deformation potential cou- 
pling, respectively. We have used 7 = 11.15 eVA 3 and 
a = for the Dresselhaus and Rashba SO coupling 
strengths, respectively. Apart from the smallest and 
largest energy gaps, where the relaxation rate is very 
small, we obtain a very good agreement with the exper- 
imental measurements. We find that interchanging the 
values between the two couplings has only a minor effect 
on the transition rate; in fact, this is true for any com- 
bination with \/ a 2 + {"f(p 2 )) 2 = 0.58 meVA. The exper- 
iment of Meunier et ali^ gives a measure of this quan- 
tity but cannot resolve the relative magnitude of the two 
terms. Measurements of angular resolved phonons 3 ^ or 
external fields controlling the Rashba term are needed to 
address the comparison between theory and experiment 
at a more detailed level. 

In general, the coupling to the phonon bath is strongest 
when the QD size matches with the phonon wave- 
length^^ ' 15 ' 27 The piezoelectric coupling (dashed line) is 
found to dominate the relaxation rate in Fig. [2] at small 
energy splittings (AE < 0.5 meV), whereas the defor- 
mation potential coupling dominates at larger AE. This 
is expected due to the V 'AE and y/T/AE dependence of 
these couplings, respectively. In the former case, the cou- 
pling occurs through slower transverse phonons, which 
yields a peak at smaller AE. 

Next we examine the effect of the SO coupling strength 
on the relaxation rate in more detail. Figure [3] shows the 
scaled relaxation rates for four different sets of a and 7 
in the units of meVA and eVA 3 , respectively. The black 
solid line corresponds to Fig. [2] The other curves are 
divided by the actual a 2 + (l(p 2 z )) 2 and multiplied by the 
square of the reference value, 7 = 11.15 eV A 3 . Overall, 
no significant qualitative changes in the rate are obtained 
within a realistic parameter range. In particular, the 
peak position is insensitive to both 7 and a, and thus we 
cannot achieve a better agreement with the experimental 
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FIG. 2: Calculated total triplet-singlet relaxation rate when 
the magnetic field tilting angle is 9 = 50° (solid line) in com- 
parison with the experimental data in Ref.— (circles with er- 
ror bars). The other parameters are the same as in Fig. [1] 
The solid line shows the total rate. The dashed and dash- 
dotted lines show the contribution from the piezoelectric and 
deformation potential coupling, respectively. The spin-orbit 
coupling parameters are 7 = 11.15 eVA 3 and a = for the 
Dresselhaus and Rashba terms, respectively. 



rate in Fig. [5] (circles) by tuning either the Dresselhaus 
or Rashba SO coupling, or both. Furthermore, the fact 
that all curves are very close to each other indicates that 
the transition rate scales to a good agreement as 

r , , [« 2 + (7(^)) 2 ] w v ncn 
r W7 ~ 2 1 7 / 2U2 r a o,7o) ■ (19) 
+ (70 (Pt)) 2 

The scaling formula can be derived from an interaction 
Hamiltonian Vj ~ aVj + 7V} 2 in the first-order pertur- 
bation limit when the interfering paths can be neglected, 
i.e., 

</|V/|i)(/|V?|i)«0. (20) 

From Fig. [3] we see that the scaling for fixed B,<j> is 
valid. This suggests that the spin relaxation through 
the Rashba or Dresselhaus coupling in fact follows non- 
intcrfcring paths within our basis states. 

Generally, however, the SO coupling is anisotropic and 
depends on the orientation of the magnetic field with 
respect to the crystal axis^ Consequently, the relaxation 
rate is anisotropic as well.— Hence, in view of a tilted 
magnetic field it is difficult to precisely assort the Rashba 
and Dresselhaus contributions in our QD device. 

In Fig. [4] we show how the relaxation rate for fixed 
magnetic field strength (here B = 2 T) depends on the 
field direction with any combination of 9 and (f>. As 
pointed out above, the relaxation rate is strongly depen- 
dent on the energy splitting and qualitatively follows a 
form AE n exp (— cAE 2 ) with positive constants n and c. 
At B = 2 T the rate is small for both 6 ~ and 6 ~ 90°, 
corresponding to small-Ai? and large-Ai? limits in the 
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FIG. 3: (Color online) Scaled relaxation rates according to 
Eq. (|19[) for four different sets of spin-orbit coupling strengths. 
a and 7 are in units of meVA and eV A 3 , respectively. 




FIG. 4: (Color online) Dependence of the relaxation rate on 
the magnetic- field angles 9 and <f>. The field strength is fixed 
to B — 2 T. The other parameters are the same as in Fig. [T] 

qualitative formula, respectively. At intermediate val- 
ues of 9 we find an area of higher relaxation rates that 
increases as a function of <f>. In this region of increased 
rates the energy splitting gives phonon wavelengths com- 
parable with QD size. Detailed analysis of the angular 
dependence is not the scope of this work, but we finally 
point out that these results could be used as a guideline 
in the analysis of forthcoming experiments where the di- 
rection of the magnetic field can be varied; in principle 
the angular properties could be used indirectly to obtain 
precise structural information of the QD device. 

C. Effects of dimensionality 

Finally we turn our attention to the effects of the di- 
mensionality on the singlet-triplet energy splitting and 



6 



on the spin relaxation rate. We compare our 3D scheme 
as described in the beginning of Sec. [II] to a bare 2D ap- 
proach, where only the in-plane confinement and 2D ba- 
sis functions are applied and the 2D approximation for 
the electron-electron interaction is used (corresponding 
to uj z 3> ujq. This leads to relatively stronger electronic 
repulsion in the ground state than in the excited states, 
which in turn yields a smaller energy splitting. 

The results of the comparison are summarized in 
Fig. El First we notice that the tilted magnetic field has a 
very different effect on the energy states in 3D and 2D (a) . 
In 2D the applied magnetic field closes the energy gap al- 
ready at B ~ 2 T, whereas in 3D the closing occurs at 
B ~ 3 T. Therefore the 2D model cannot yield a reason- 
able agreement with the experimental splitting as shown 
in Fig. EJb) . It should be noted that if the confinement 
strength ujq in the 2D model is slightly increased to yield 
the correct splitting at B = 0, the magnetic-field depen- 
dence is still wrong, and the rate becomes much worse 
than the 3D result in comparison with the experiment. 
Furthermore, as shown in Fig. [SJc) the 2D model is un- 
able to reproduce the experimental spin relaxation rate. 
In contrast, the 3D model has a very good agreement 
with the experiment as discussed above within Fig. [2] 



IV. SUMMARY AND DISCUSSION 

To summarize, we have proceeded towards quantita- 
tive modeling of spin relaxation in semiconductor quan- 
tum dots containing two electrons. Our approach covers 
(i) a three-dimensional description of the quantum dot 
device, (ii) numerically exact treatment of the electron- 
electron interaction, (iii) a tilted magnetic field, (iv) the 
ellipticity of the quantum dot, and (v) both the Rashba 
and the Dresselhaus spin-orbit coupling, and the latter 
with both linear and cubic terms. We have attributed the 
observed nonlinear behavior of the singlet-triplet energy 
gap at B < 0.5 T to the ellipticity of the quantum dot. 
Then we have found a good agreement between theory 
and experiment in the triplet-singlet energy splitting and 
in the relaxation rate, although the optimal tilting angle 
of the magnetic field was found to be considerably smaller 
than the experimental value. Finally we have explicitly 
shown that the three-dimensional model is essential to 
obtain a reasonable agreement with the experiment in 
both the singlet-triplet energy splitting and the spin re- 
laxation rate. 

In view of our results supplied with aspects (i)-(v) 
listed above, the origins behind the remaining quanti- 
tative deviations between experiment and theory - espe- 
cially regarding the different angle of the magnetic field 
- is unclear. However, as the most probable scenario we 
may suggest that the actual confining potential deviates 
from harmonic, e.g., by being more strongly confined at 
the edge of the quantum dot. This can be induced by the 
gate geometry of the lateral device and/or the nonuni- 
formity of the confinement in the vertical (z) direction. 
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FIG. 5: (Color online) Comparison between two-dimensional 
and three-dimensional calculations, in both cases with a tilted 
magnetic field of = 50°. (a) Energy levels of the lowest 
singlet and triplet states, (b) Singlet-triplet energy splitting, 
(c) Relaxation rate. Parameters are the same as in Fig. [2] 



These conditions can increase the high sensitivity of the 
spectrum to the tilted magnetic field (Fig. [I} even fur- 
ther. In this respect, we hope that the present study mo- 
tivates more efforts in this direction, both theoretically 
and experimentally. 
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